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We use the coupled cluster expansion (exp(S') method) to generate the complete ground 
state correlations due to the NN interaction. Part of this procedure is the calculation of 
the two-body G matrix inside the nucleus in which it is being used. This formalism is being 
applied to ^^O in a configuration space of 50 huj. The resulting ground state wave function 
is used to calculate the binding energy and one- and two-body densities for the ground state 
of 1^0. 
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I. INTRODUCTION 

In the last thirty years electron scattering from nuclei 
has provided a wealth of information mapping out nu- 
clear ground state charge densities [1], providing precise 
transition charge and current densities for the excitation 
of single particle states [2] and for collective states [3]. 
The measurement of ground state magnetization densi- 
ties and the excitation of high multipolarity magnetic 
excitations, or the single particle knockout reaction to 
discrete states all have in some way supported the mean- 
field approach as the lowest order in the description of 
nuclear structure. 

The confirmation of the mean-field approach, however, 
were more qualitative in nature than quantitative. The 
form factors e.g. for the excitation of the high spin sin- 
gle particle states in ^°^Pb [4], were described extremely 
well in shape by the mean-field wave functions, however, 
the predicted strength was too big by a factor of two. 
The knockout reactions again were in good agreement 
with the shapes predicted by the mean-field wave func- 
tions but the strength was off by again roughly a factor of 
two [5] . The general conclusion was that correlations are 
important in the calculation of observables. They do not 
change the shape of the wave functions but they mod- 
ify the strength due to deoccupation of orbits below the 
Fermi surface and partial occupation of the orbits above 
the Fermi surface [6]. This was confirmed by (e,e'p) ex- 
periments in which particles from orbits above the Fermi 
level were knocked out [5]. Thus, to do justice to the ac- 
curacy of the electromagnetic probe, we can no longer be 
satisfied with the mean-field approach but have to take 
into account the correlations largely due to the hard re- 
pulsive core of the nucleon-nucleon interaction. While we 
put particular emphasis on the effects of correlations on 
the results of electron scattering experiments this is less 



of an issue in this paper dealing with the ground state, 
but will be more clearly pointed out in subsequent papers 
dealing with excited states of these nuclei. 

There are different ways to account for correlations. 
One way is to introduce correlation functions in the many 
body wave function in real space. This has been quite 
successful for small nuclei [7] and has resulted in rea- 
sonable descriptions of the ^^O nucleus [8]. A different 
approach is to add in configuration space to the uncor- 
related ground state multi-particle multi-hole configu- 
rations. Both approaches can be related to each other 
even though it seems harder using correlation functions 
to satisfy the Pauli principle. In configuration space this 
appears to be more transparent. In our treatment we 
follow closely the formulation of the Bochum group [9]. 
However, we solve the equations entirely in configuration 
space. Furthermore, we truncate in different ways where 
the significance of terms becomes more transparent 



II. THE UNCORRELATED GROUND STATE 

We assume an orthonormal set of single particle wave 
functions exists. These wave functions are solutions to 
the single particle Hamiltonian given by the Schrodinger 
equation in the mean-field. This mean-field will be gen- 
erated in an iterative procedure. The single particle wave 
functions are expanded in a set of orthonormal functions 
such as harmonic oscillator functions, Bessel functions 
or others. Thus each orbit is represented by a set of 
expansion coefficients and a single particle energy. An 
uncorrelated ground state can be constructed as a single 
Slater determinant which includes all the occupied orbits 
and is written as |0), our new vacuum state. 



III. VARIATIONAL CALCULATIONS FOR THE 
CORRELATED GROUND STATE 
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We use the formalism of second quantisation. As such 
we follow the notation of the textbook by deShalit and 
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Feshbach [10]. In the second quantisation representation 
we write the Hamiltonian of a closed shell nucleus as 



H = Tfcife2 4i + 



, ^lfe2fc3fe4 4i 



afe4 afeg 



(3.1) 



where the Einstein summation convention is understood 
and Vk-^k2k3ki are fully anti-symmetric matrix elements of 
the potential. The orbits completely occupied in a non- 
correlated ground state are called "hole orbits", those 
not occupied are called "particle orbits". Creation and 
destruction operators in these orbits are denoted a^, aj^, 
and ap, a/j, respectively. Thus, we have 



al |0) 







|0) = 0. 



(3.2) 



The correlated ground state is assumed to be of the 
form 



|6) = e^' |0) = (1 + Ft) |0) . 



(3.3) 



Here is the operator that describes the correlation. It 
is written as 



(3.4) 



where 



ph 

^2 ~ $Z ^Pihi,P2h2 ^pi ^p2 ^h2 ^hi ) 

Plhi,p2h2 



(3!) 



(3.5) 



or correspondingly F^. The purpose of the variational 
calculation is to determine all the coefficients Z. 

We will further assume that there exists at least one 
set of wave functions, satisfying the condition that Si =0. 
Further down we will discuss how to construct such a set. 
This set is the set of single particle mean-field wave func- 
tions. With that basis the lowest order of correlations 
are the 2p2h correlations. Also, S2 and F2 are the same 
as well as S3 and F3, whereas F4 = S4 -|- S2/2, etc. 

A variation S\0) orthogonal to the correlated ground 
state can be constructed from any operator Cj, repre- 
senting any npn/i-excitation as 

^|6) = e-^cte-s* |6) = e-^Ci |0) . (3.6) 



We have 



<5|6) = (0| ct |0) = 0. 



(3.7) 



The variational principle requires that the Hamiltonian 
between the ground state and such a variation vanishes. 
Thus, we have 



(0| H (5|0) = (0| H e"^ |0) = . 



(3.8) 



Choosing C], to be Iplh, or 2p2h, ApAh results in a 
set of nonlinear coupled equations that are written out 
in more detail in the next section. We might say that 
e^He~ ^ represents the effective Hamiltonian as Eq. (3.8) 
represents the Hartree-Fock condition for the uncorre- 
lated ground state with this effective Hamiltonian. 



A. The effective one-body Hamiltonian 

To simplify the accounting of number of p/i-excitations 
we use indices for the Hamiltonian. Explicitly we define 



V20 
Vio 

Voo 



-^VpiP2hih2^pi^P2^h2^hi — , 



~ 2^PiP2hiP3^i^P2^3^hi 

+ 2^Pih2hih3^li^h2^'^3^hi = V| 

V^i /i2?iiP2 ajj^ ap2 a/ij 

1 _ II 

^ ''P1P2P3P4 ^pi ^P2 ^Pi ^P3 



- Vt 



01 ' 



1, 



+ ^Vhih2h3h4^ln^h2^h4^h3 

Tkik2 a^j^a/cj . 



To 



(3.9) 



We should take note that any S„ commutes with Vo2- 
We introduce the mean-field Hamiltonian Hq as 



Ho = To + Uo 



(3.10) 



where To is the kinetic energy operator and Uo is the 

mean-field to be specified below. We will assume that the 
orbits are eigenfunctions of this mean-field Hamiltonian 
Ho with 



[Ho, at] 



e^at 



[Ho, Ah] = -ChSih ■ (3.11) 



We use Eq. (3.8) with n — \, C\ = ata/j, where we 
write ata/j |0) = \lplh) to obtain the equation establish- 
ing Si 

= (0||to + Vol + [Si, To] + [Si, Voo] + [S2, Vio] + 

Si, [S2,V2o] 



[S3,V2o] + 2[Si, [Si,Vio] 



+ 



Si, 



Si, [Si,V2o] 



\lplh) . 



(3.12) 



There are similar equations that determine S2,S3,.... 
While these equations hold in any basis, there is one basis 
of particular convenience. This is the maximum overlap 
basis in which Si vanishes. Equation (3.12) results in 
the solution Si = if the terms that do not contain 
Si vanish. The mean- field basis is determined by the 
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condition of the vanishing of Si and in the mean- field 
basis we must have 



+ V, 
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0= (0| |(Ho-Uo) 

+ [S2,Vio] + [S3,V2o] I \lplh). (3.13) 

Using Eq. (3.11), wc can show that the expectation value 
(0|Ho|lpl/i) vanishes. Therefore, Eq. (3.13) becomes 

(0|Uo|lpl/i) = (Oljvoi + [S2, Vio] + [S3, V20] ||lpl/i) . 

(3.14) 

Thus those terms establish the elements in the one- 
body Hamiltonian matrix that connect p and h orbits. 
The equations establishing the higher order correlations 
in the mean-field basis are 

= (0||vo2 + [S2, Voo] + [S2,Ho] - [S2,Uo] 

+ [S3,Vio] + [S4,V2o] +i[S2,[S2,V2o]]||2p2/l), 

(3.15) 

= (0||[S2,Voi] + [S3, Voo] + [S3, Ho] - [S3,Uo] 

+ [S4, Vio] + [S5, V20] + ^ [S2, [S2, Vio] 

+ [S3, [S2,V2o]]||3p3/i), (3.16) 

0= (0||[S3,Voi] + [S4,Voo] + [S4,Ho] - [S4,Uo] 

+ [S5,Vio] + [S6,V2o] + ^[S2, [S2,Voo] 



+ 



1 



S3, [S2, Vio] + S4, [S2, V20 

1 



S3, [S3, V20] 



+ 



'2, 



S2, [S2,V2o] 



|4p4/i) . 
(3.17) 



At this point we will assume that the orbits arc cigen- 
functions to the single particle Hamiltonian Hq. This 
allows us to solve these equations as 



|S2|2p2ft) 



V02+ [S2,Voo] - [S2,Uo] 



+ [S3, Vio] + [S4,V2o] 

1 

+ 2 



S2,[S2,V2o] \ — \2p2h), (3.18) 



S3|3p3/i) = -(0|<^ [S2,Voi] + [S3, Voo] - [S3,Uo] 



+ [S4,Vlo] + [S5,V2o] + 2 [S2, [S2,Vio] 



(3.19) 



S3, [S2,V2o]l \^J^p3h) 



A similar equation allows us to isolate S4 using 
Eq. (3.17). 

If we restrict ourselves to a situation in which only 

2p2h correlations arc present, S2 and F2 become identi- 
cal. Then, Eq. (3.18) can be written in its adjoint form 
as 



vfI). 



(3.20) 



This equation is essentially the Bcthc-Goldstonc equation 
for the G-matrix which is normally written as [12] 



Gt(a;) = V 



1 



Ho 



(3.21) 



Multiplying this equation from the left with l/(w — Ho), 
using F2 = [w — Ho)]""^ G^(a;) and setting w = 0, results 
in Eq. (3.20) for Thus, one of the essential parts 
of the coupled cluster expansion is that we calculate the 
ground state two-body G-matrix inside the nucleus where 
it is to be applied. The additional terms in our coupled 
equations indicate the need for appropriate corrections 
in the G-matrix for the presence of 3p3h, ApAh, etc. cor- 
relations. In most cases the G-matrix is calculated in 
nuclear matter and then applied to finite nuclei. How- 
ever, in nuclear matter the correlation function extends 
up to infinity and thus cannot be applied to a finite nu- 
cleus. 

Wc estimate that in our basis there arc about 10^ 
2p2h configurations and about 10^° 3p3/i configurations. 
While the number of 2p2h configurations is quite acces- 
sible, the number of ZpZh configurations is prohibitively 
large, and we cannot store all these numbers. Thus we 
have to implicitly correct for the presence of these cor- 
relations. Wc do this by inserting the solutions for S„ 
with n>3 back into the equations and thereby obtaining 
a perturbation expansion in l/Eph- We write this equa- 
tion out up to second order for the Eq. (3.14) establishing 
the mean-field 



(0|Uo|lpl/i) 



|Voi|lpl/^) + (0|[S2,Vio]|lpl/i) 



-(0| [S2,Voi|,— V20 



:(0| 



+ 



S2,[S2,Vio] 

1 



\lplh) 

1 



Ho 



V. 



20 



\lplh) 



[S2,VoiJ,:z^Voo 
■Ho 



1 

H^ 



V20 



\lplh) . (3.22) 



This equation establishes the matrix elements of the 
single particle Hamiltonian Hq between particle and hole 
orbits. The matrix elements between hole and hole orbits 
or between particle and particle orbits are not defined. 
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and any definition may be chosen. As long as Uo is ex- 
plicitly kept on the right hand side of Eqs. (3.18, 3.20) 
the explicit choice is merely a question of how fast the re- 
sulting series will converge. However, a reasonable choice 
appears to be that form that we obtain if we replace in 
the matrix elements obtained in Eq. (3.22) the hole orbit 
with a particle orbit in order to get the matrix elements 
between particle and particle orbits and we change the 
particle orbit into a hole orbit in order to got the matrix 
elements between hole and hole orbits. Reference [13] 
gives a detailed account of the contributions included in 
our mean- field as given by Eq. (3.22). Our choice for the 
other matrix elements corresponds simply in turning the 
hole line into a particle line or vice versa. 

The mean-field orbits are the eigenvectors of this ma- 
trix, and the eigenvalues are the single particle energies. 
This procedure now fully defines the mean-field used here 
even though its definition is not unique. 

To obtain values for the amplitudes of S2 we start with 
Eq. (3.18) and replace S3 by Eq. (3.19) and correspond- 
ingly S4. We continue with these replacements, resulting 
again in an expansion in 1/ Eph for S2 . 

The calculation of the correlations breaks down into 
two steps: In step one the two-body correlations arc 
computed according to Eq. (3.18) where the curly 
brackets give the various contributions to the effective 
{pihi\Vph\h2P2) matrix element. In step two the single 
particle energy tensor is calculated from the kinetic en- 
ergy, the direct potential energy, and the correlation en- 
ergy. The eigenvalues and the eigenvectors of the energy 
tensor give the new single particle energies and the now 
single particle functions. The two steps are iterated un- 
til a self-consistent solution is obtained, i.e. the energy 
tensor is diagonal. 

As wc use an infinite expansion we must truncate that 
expansion. We find that terms connecting to effective 
one-body terms give the major contributions, two-body 
terms are less important. Generally, we have included 
all terms that arc written as product of three two-body 
operators. Of those terms written as product of four 
two-body operators we have included those terms where 
at least two operators connect to a one-body operator. 
This results in all the quenching terms for the products 
of two two-body operators. We have also included those 
terms that look like the terms of the products of three 
two-body operators with a renormalized interaction. We 
found that these corrections are equally significant for 
those terms arising from S4 and those arising from S3. 
Left out are mostly those terms that required several 
angular momentum recouplings. Contributions from S5 
show up only if one considers products of more than four 
two-body operators. In the Ref. [13] we have listed ex- 
plicitly the relations and approximations that we have 
used in this calculation. 



IV. NUMERICAL CALCULATIONS FOR ^^O 

We have solved the main Eq. (3.22) that determines the 
2p2/i-amplitudes and thus essentially the ground state 
G-matrix for ^^O in a space of 50 htu with a harmonic 
oscillator length parameter b=0.8 fm, excluding those or- 
bits with I > 13. Corrections for 3p3h correlations were 
included in a reduced space of 30 Hlo and / < 6, while 
correlations due to 4p4/i correlations were included in 
the full space. We used the Argonne vis potential [11] to 
generate the matrix elements [14]. For comparison, re- 
sults for the Argonne and Vu potentials arc reported. 
Here, the Argonne vs potential is the reprojection of the 
Argonne V14 potential in the sense of reference [7] . 

The Hamiltonian is given in the center-of-mass as 



A ^ A 



(4.1) 



where Tcm = P^/2M is the kinetic energy operator of 
the center-of-mass, and M is the total mass. This repre- 
sents the energy in the center-of-mass frame. It can be 
rewritten as 



Hi, 



1 

~ A 

A 

+ E 

i<j=l 



1 



A 

E 

1=1 



Pi ■ Pj 

Am 



(4.2) 



The last term is treated as part of the internal potential, 
and the antisymmetric matrix elements of the internal 
potential include this center-of-mass correction term. 



A. Ground state expectation values for arbitrsiry 
operators 

Ground state expectation values can be evaluated by 
introducing the operator S\ as described e.g. in the 
review by Bishop [12]. We take the ground state as 



|0) = e^' |0) . 



(4.3) 



We have defined {Cjj} as the complete set of Iplh, 2p2h, 
ApAh excitations. The normalized expectation value 
a of any operator A, a = (A), can be worked out as 



(0| A e^^ |0) _ (0| A e'^ e^' |0) 



a 



(0|0) 



(0|0) 



(4.4) 



By inserting the unity operator in the {Cj^ |0)} basis, we 
obtain 
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a = (0| A e"^ |0) 

+ E<OI«=A«-=Cno><5^^i£fjf^. (4.5) 

The expectation value on the right is by definition c„, 
the expectation value of C„. Thus we can define the new 
operator 



(4.6) 



n=l 



With this, the expectation value for any operator can be 
expressed as 

a = (0| A e-^ (l + St) |0) . (4.7) 

The operators can be obtained by solving Eq. (4.6) in 
an iterative fashion. Explicitly we write in the same 
form as Eq. (3.5), but use the amplitudes S instead of Z 
to distinguish it from S^. 



B. Ground state binding energy 



equations we have ignored some of the couplings back 
into S2. As a result Eq.(4.11) is no longer exact and the 
feature of being an upper bound to the true energy is lost. 
We believe that the errors due to the truncations are rea- 
sonably small. However, more experience is needed with 
some of the terms left out in order to reduce the uncer- 
tainties in the result. Table I shows the resulting binding 
energy for the Argonne vg, fi4 and ?;i8 potentials. 

Our Hilbert space is controlled by two cut-off param- 
eters Imax and Umax ) which are related to the size of the 
Hilbert space as N^ax = ^^max + Imax- We have stud- 
ied the dependence of the binding energy on these two 
parameters. Figure 1 displays the Imax dependence and 
shows reasonable convergence for Imax — 

11. This de- 
pendence was mapped out with nmax='^5. We have also 
checked the sensitivity of the binding energy with respect 
to rimax with a fixed Imax- In this case, we find that for 
our choice of a 30 Hlo 3p3h space, the binding energy is 
independent of rimax for an rimax value between 22 and 
25. As we go to larger Umax values however, we start 
seeing the effect of having a smaller 3p3h space and the 
approximation starts breaking down. The dependence 
of the binding energy on the Umax cut-off is depicted in 
Fig. 2 for lmax='i-'i- 



We first apply the above formalism to the ground state 
binding energy. The expectation value of the internal 
Hamiltonian can be written as 

e = (0| H,„i e-^ |0) . (4.8) 

Because of the Hartree-Fock (HF) condition expressed in 
Eq. (3.8) the terms involving vanish and we get 

e = (0| H e"^ |0) . (4.9) 

Assuming that H is at most a two-body operator and 
taking into account that Si vanishes, we write this as 

{E) = e = (0| H |0) + (0| S2 V20 |0) . (4.10) 

When we evaluate the expectation values of the opera- 
tors in the above equation, we have to consider that the 
hole orbits are not diagonal with respect to any of these 
operators. Also, as we have used the HF condition this 
expression does not give an upper limit of the ground 
state energy unless we are exactly at the minimum. In 
terms of matrix elements the energy can be written as 

h\h2 hih2 

+ ^ ^ ] ^PiP2,hih2 ypiP2,hih2 - (4-11) 
PiP2hih2 

The use of Eq. (4.11) for the energy implies that all the 

correlations arc present and satisfy the Eqs. (3.8). In 
that case the resulting energy could be taken as an up- 
per bound to the true energy. However, in solving these 



C. Ground state one-body density 

Next, we apply this procedure to the ground state 
one-body density. The expressions necessarily look more 
complex as we cannot apply the HF condition to simplify 
the expressions. 

By definition, the ground state one-body density is in- 
troduced as 

A 

p(r) = ^(<5(r-rl)). (4.12) 

Since we are dealing with a spherically symmetric system, 
we shall integrate out the angular degrees of freedom of 
the system. Then, we write the density operator in the 
second quantisation representation as 

p-Pir) = E Pa/3(r)aU/3- (4-13) 

a/3 

Here we use pa/3 (r) = Ra (r) Rf^ (r) to denote the radial 
part of the expectation value (a | 5{f—r') | Thus, we 
can write the ground state density as 

p{r) = E Rair) Rp{r) , (4.14) 

a/3 

where 

da,b = ( a/3 ) (4.15) 

is the density matrix. The density matrix is a real sym- 
metric matrix with positive definite eigenvalues. We can 
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make a basis transformation such that the density matrix 
becomes diagonal. This basis represents the "natural" 
orbits. In this basis the density becomes 

P{r) = E [RT\r)Y . (4.16) 

a 

Here v™* represents the occupation probability of these 
natural orbits. This is the only basis in which occupation 
probabilities have a moaning. 

To calculate the one-body density matrix we use 
Eq. (4.15) with Eq. (4.7) 

d^p = (0| ata/3 |0) + (0| [S2, a^a/j] S| |0) 
+ (0|aia^Sl |0) + (0| [S2, ata^] sj 1 0) 
+ (0| [S3, ala^] S| |0) + (0| [S3, a^a^] Sj 1 0) 

+ i (0| [S2, [S2, a^a^j]] Sj |0) + ■•• . (4.17) 

In evaluating terms of equal order we have to consider S2 
to be of order 1/Eph, S3 to be of order l/E^/^, S4 to be 
of order l/E^f^, etc. To account for S3 and S4 wc have 
to use the expansions of Eq. (3.18). As Si vanishes, Si 
is of the same order as S3. Otherwise S„ is of the same 
order as S„. 

The one-body density calculated this way is not the 
density in the center-of-mass frame of the nucleus as the 
wave function represents a nucleus that has a residual 
center-of-mass motion. The form factor measured in ex- 
periments is the form factor in the center-of-mass frame 
which we can write as 

f{q) = ^ ( e'«-(*^'-'^-'") ) = ( e'^-^^ ) { e"*?-^-- ) . 
j 

(4.18) 

i.e. the form factor calculated from the one-body den- 
sity factorizes into the form factor in the center-of-mass 
frame and the form factor of the center-of-mass motion 
Fcm{q) = (e-*?-«=-). 

It has been shown that for a pure uncorrelated har- 
monic oscillator wave function this form factor can be 
written as Fcmiq) = exp{—b^q^ /AA) where b is the har- 
monic oscillator length parameter. There are two rea- 
sons why this docs not apply here. First, wc have self- 
consistent mean-field wave functions and not harmonic 
oscillator wave functions and second, we have not a sin- 
gle Slater determinant but a sum over many due to the 
correlations. For these reasons we chose to evaluate 
the form factor directly using Eq. (4.18). By writing 
Rem = ^ X^i^i) the Eq. (4.18) is expanded into n-body 
terms. This expansion was checked for the case of a sin- 
gle Slater determinant of harmonic oscillator functions 
against the exact result. It was found that this expan- 
sion is quite satisfactory if all terms up to three-body 
terms are kept [15]. 



Fig. 3 shows the calculated charge density after folding 
the proton point density with the charge density of the 
proton and folding the neutron point density with the 
charge density of the neutron. As this expansion is accu- 
rate up to terms of order g^, it encompasses the result of 
the rms-radius. At the present time no corrections due 
to the meson-exchange charge density are taken into ac- 
count. The resulting charge radii are shown in Table I 
for the Argonne vs, Wi4 and fig potentials, respectively, 
and are reasonably close to the experimental one. 

In the calculation of the natural orbits we also gener- 
ate the occupation probabilities for the orbits above the 
Fermi level. For the Argonne vs, wi4 and vis potentials 
the occupation probabilities of the 1^5/2 and the 2si/2 
proton orbits are summarized in Table I and appear to 
be consistent with the experiment [5], which establishes 
a lower limit for these values. 



D. Ground state two-body density 

A direct presentation of the short range correlation can 

be seen in the ground state two-body density. We start 
with the ground state two-body density definition 

p(ri, rs) = 5Z <0 I -^(^i - "^(^2 - fj) 1 0) . (4.19) 

In the second quantization representation the two-body 
density operator can be written as 

^''^(ri, f2) = E (a/3 I p(ri, ra) | j6} a^ a^ a^ a^ . 

a/375 

(4.20) 

Using the completeness relationship of the spherical har- 
monics we can evaluate the matrix element Pai3-yS 

{aP\p{fi,f2)\jS) 

2'm2 I .?7^7/ 

/imi 

X XI R(3ir2) Rsir2) Yi^rn2{f2) {jfjmff I Yil^^ I jsms) . 

I2m2 

(4.21) 

In order to be consistent with the phase convention of 

the two-body potential matrix elements, wc couple the 
two-body density matrix elements using the ph. angular 
momentum coupling conventions [14] . Then, the angular 
momentum coupled density is 

((a7)A|p^''(ri,r2)l(^^)A) 

= pU^i) pU^2) lA;(ri) n^(f2) . (4.22) 

Here we have introduced the one-body multipole density 
p^^(r) which is 
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= ^(2i„ + l)(2j> + l) 

0,1/2 j> - 1/2 I AO) i?„(r) ii;3(r) (4.23) 

it la + IjS + ^ is oven, and zero otherwise. 

For a spherically symmetric (spin=0) nucleus it is more 
relevant to calculate r2, 6*12) as due to the spherical 
symmetry the two-body density is dependent on the di- 
rection of fi relative to alone. Thus, we can perform 
an average over the directions of r2- This translates into 
carrying out the sum over the /i component of the angular 
momentum A. We obtain the result 

Pa0^s{ri,r2,Oi2) = pt-yin) Ps/3{r2) Px{cosei2) (4.24) 

We now discuss the most dominant contributions. We 
again apply Eq. (4.7) to evaluate the two-body density 
matrix, ( aj^ a^ a^, ). With this, we get the ground 
state two-body density as 

p{ri,r2,0i2) 

= {o\p"/\o) + {o\pTsI\o) + (o|S2prio) 
+ (0|[S2, pT] sl\o) + (oip^^slio) 

+ i(0|[S2, [S2, p^^]]St|0) + .... (4.25) 

(a) Contributions from bare ground state: 

p{ri,r2,9i2) 

= E + 1) RlM) (2i/^. + 1) «L(^2) (4.26) 

hih2 

This term has the structure phf{l)phf{'2)- Here phf is the 
bare one-body density. We can include all higher order 
one-body density terms by adding 

Ap{ri,r2,0i2) = Ap{ri)phf{r2)+phf{ri)Ap{r2) , 

(4.27) 

where Ap(r) is the difference between the full one-body 
density minus the bare one-body density. In turn, in 
higher order terms we have to exclude those terms where 
in the operator Paj3-yS ct connects with 7 or /3 connects 
with 6. However, we still have to consider the exchange 
terms where a connects with S and j3 connects with 7. 

The exchange term is due to the Pauli correlations and 
results in 

Ap{ri,r2,0i2) 

= - E pUM) Plh.ir2) Pl{cose^2) ■ (4.28) 

hih2 I 

(b) Contributions linear in the 2p2/i-amplitudes. For 

the following terms we assume a summation over all or- 
bits appearing twice in the expression. The leading term 
that arises from the correlations is 



Ap(ri,r2,^12) = E iSp^h^p2h2 + Z^^hrp^h^) 
A 

xplpA^i)pU.(^2)Px{cose,2) . (4.29) 

In our application we have included all terms up to 
second order. Aside from the terms listed here, the next 
most important correction arises from the second order 
term depending on S2 and S^- 

The two-body density represents the probability of 
finding one nucleon at fi and one nucleon at 7%. We can 
divide this by the probability of finding the first nucleon 
at fi. The remaining density represents the probability 
of finding a second nucleon at r2 if the first nucleon is 
at fi. If both nuclcons arc protons, this density is nor- 
malised to a total integral of (Z-1). Figs. 4 to 9 show 
these densities as a function of r2 for various positions 
of ri. We have made no attempt to correct these for 
the residual center-of-mass motion of the nucleus. The 
densities show the effects of the short range repulsion: 
they exhibit a deep hole where the first nucleon is lo- 
cated. The two-body densities also show that for large 
distances the long-range aspect of the ground state nu- 
clear correlations, usually thought to be related with the 
surface deformation modes, has a significant contribu- 
tion: when the first nucleon is located closer to the nu- 
clear surface, we observe an enhancement of the density 
at the symmetrically-opposite position. The picture of a 
two-body density obtained as the revolution of the spher- 
ically symmetric one-body density, with a Gaussian-like 
distribution centered at the location of the first nucleon 
scooped out of it, is definitely insuSicient. 

V. CONCLUSIONS 

Wc use the coupled cluster expansion (exp(S') method) 
to solve the many-body Schrodinger equation in configu- 
ration space. While the coupled cluster expansion is ex- 
act if carried out to all orders, the present results are ob- 
tained with truncations. As we have retained only terms 
up to second order in 1/e the contributions from S5 have 
no effect on 82- In turn, the effects of any higher n cor- 
relations are only of the order as the higher order terms 
left out anyway. Thus, in essence we have replaced the 
truncations in S„ by the more relevant truncation in 1/e. 

We have shown that it is possible to choose large 
enough configuration spaces for the complete and self- 
consistent calculation of the ground state correlations in- 
side a finite nucleus. This calculation makes no artificial 
separation between "short range" and "long range" cor- 
relations. In fact, the two-body density shows that the 
correlation function in the surface region of the nucleus 
has strong contributions from the surface deformation 
modes. It is largely these modes that cause the strong 
deoccupation of orbits close to the Fermi surface. The 
occupation of the 2si/2 and the lc?5/2 proton orbits cal- 
culated is consistent with the values suggested by the 
experiments. 
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Our efforts arc currently directed in two directions. 
First, we intend to apply the procedure described in this 
paper to the particular case of a more realistic interac- 
tion, namely the Argonne Vis potential [11] together with 
a phenomenological three- nucleon interaction [16]. This 
should not only result in a better description of the ^^O 
observables, but a breakdown of the contributions to the 
binding from the two- and the three-body interactions 
could be inferred also. Secondly, we shall use the equa- 
tion of motion technique to calculate excited states of the 
^^O nucleus, and the ground state and excited states of 
the neighboring nuclei (e.g. ^^N). We hope to be able to 
present our findings in the near future. 
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TABLE I. Energy expectation values, charge radii, and 
proton orbits occupation probabilities, for the Argonne vs, 
1714 and Vis potentials, respectively. 



Potential 



vU 
expt. 



E 

[MeV/nucleon] 



rms 
[fm] 



7.0 
6.1 
5.9 

8.0 



2.81 
2.86 
2.81 

2.73 
± 0.025 



Id: 



■5/2 



3.68 
3.33 
2.58 

2.17 
± 0.12 



2Si/2 

[%] 



4.09 
3.99 
2.75 

1.78 
± 0.36 



-4.0 
-4.4 
-4.8 
E -5.2 
-5.6 
-6.0 
-6.4 
-6.8 



10 
max 



FIG. 1. Dependence of the binding energy on the Imax 

cut-off for nmaa; =25. 



22 23 24 ^25 



FIG. 2. Dependence of the binding energy on the Umax 



cut-off for 



=11. 




1 2 3 4 5 



r[fm] 

FIG. 3. Comparison of the "experimental" with the calcu- 
lated charge density, for the Argonne vs, vu and vis poten- 
tials, respectively. 
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FIG. 4. The p-p two-body density for the Argonne vis 
potential: the case when the first proton is located at 
xi = 0.0/m 




FIG. 5. Same as Fig. 4 with the first proton located at 
xi = 1.0/m 




FIG. 6. Same as Fig. 4 with the first proton located at 
xi = 2.0/m 




FIG. 7. The p-n two-body density for the Argonne wis po- 
tential: the case when the proton is located at xi — 0.0/m 




FIG. 8. Same as Fig. 7 with the proton located at 
xi = 1.0/m 




FIG. 9. Same as Fig. 7 with the proton located at 
xi = 2.0/m 
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